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ABSTRACT 

An unexpected space debris population has been detected in 2004 lSchildknecht et al. I 
(|2003, 2001) with t he uniq ue properties of a very high area-to-mass ratio (HAMR) 
ISchildknecht et al I (|2QQ5aD . Ever since it has been tried to investigate the dynamical 
properties of those objects further. The orbits of those objects are heavily perturbed 
by the effect of direct radiation pressure. Unknown attitude motion complicates orbit 
prediction. The area-to-mass ratio of the objects seems to be not stable over time. 
Only sparse optical data is available for those objects in drift orbits. 

The current work uses optical observations of five HAMR objects, observed 
over several years and investigates the variation of their area-to-mass ratio and 
orbital parameters. A normalized orbit determination setup has been established and 
validated with two low and two of the high ratio objects, to ensure, that comparable 
orbits over longer time spans are determined even with sparse optical data. 
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1 INTRODUCTION 

T* The Astronomical Institute of the University of Bern 
^ ' (AIUB) detected high are a-to-mass ratio (HAMR) objects 
^ ! in GE O-like orbits in 2004 ISchildknecht et al] (|200l 12004 
^ ■ l2005al ). Since then, the AIUB observes HAMR objects 
on a regular basis and keeps a small catalog of HAMR 
and other space debris objects, which are not listed in the 
USSTRATCOM catalog. The observations are performed 
with the one meter ESA Space Debris Telescope (ESASDT), 
located on Tenerife, Spain, and the one meter Zimmerwald 
Laser and Astrometry Telescope (ZIMLAT), located in 
Zimmerwald, Switzerland. Additional observations for some 
objects, which were detected by the AIUB, are provided by 
courtesy of the Keldysh Institute of Applied Mathematics, 
Moscow, via the ISON network. 

Maintaining a catalogue of HAMR objects is especially 
challenging due to the unique properties of these objects; 
the orbits are highly perturbed by direct radiation pressure. 
Regular observations on short time intervals are mandatory. 
In routine orbit determination for catalogue maintenance, 



variations in the value of the effective area-to-mass ratio 
(AMR) were detec ted, first investigations were performed 
iMusci et al7\ (l2010l ). 

For the investigations presented in this paper orbits 
are determin e d with an enhanced version of the CelMech 
tool iBeutlerl (l2005l ). The area-to-mass (AMR) value is 
determined as a scaling parameter of the direct radiation 
pressure. The acceleration due to the direct radiation 
pressure is calculated as: 
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where f is the geocentric position of the satellite, 
the geocentric position vector of the sun, AU the as- 
tronomical unit, A the effective cross section exposed 
to the radiation, m the mass of the satellite, and c the 
speed of light. C is the refiection coefficient. The direct 
radiation pressure is determined under the assumption of 
a spherically shaped object. In contrast to the calculation 
of the radiat ion pressure a cceleration by other sources 
(compare e.g. IVallado Mc Cain (200ll)), the coefficient C 
is divided by two in the formula above. A value for C has 
to be chosen, by default, 2.0 is selected in the standard 
processing. This corresponds to an assumption of full 
absorption. All AMR values presented in this paper have 
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Table 1. Internal name, eccentricity, inclination (deg), 
semi-major axis (km), area-to-mass ratio {rn^/kg) and 
apparent magnitude (mag) of the selected objects of the 
AIUB catalogue 

NAME Epoch a e i AMR Mag 



observations are allowed per set. An observation set may 
consist of more than one tracklet. But the observations 
within the sets should not be distributed over more than 
three days. 

Orbits were determined for different spacings of two 
observation sets stemming a) from one observation site only 
and b) from different sites. In the first case, the observations 
either stem from ZIMLAT or from ESASDT only. In the 
second case, not only the observations of ZIMLAT and 
ESASDT were combined but also observations of the ISON 
network, if available. When observations from different sites 
are used in orbit determination, the distribution is either 
that the first set of observations stems from one site and the 
second from another, or that there are observations from 
different sites at similar epochs used within the first and/or 
the last set of observations or a mixture of those options. 
In the figures the label ALL is applied, when observations 
of ZIMLAT (labeled ZIM), the ESASDT and of the ISON 
network are combined; the label SDT-ZIM is applied, if 
only the observations of ZIMLAT and the ESASDT are 
used. The distances between the observations and the 
ephemerides of the predicted orbits of the four objects for a 
prediction interval of 50 days after the last observation used 
for orbit determination were determined. The distances 
were averaged and a mean value and standard deviation 
was calculated. Between six and 50 single distances between 
ephemerides and observations were averaged. 

The predicted ephemeris positions are compared to 
the optical angle-only observations, which were not used 
in orbit determination. Angular distances are determined 
on the celestial sphere. The observation used for the 
comparison stem from ZIMLAT and ESASDT and serve as 
ground truth. Calibration measurements with high accuracy 
ephemerides of Global Navigation Satellite System (GNSS) 
satellites provided by International GNSS Service (IGS) 
showed an accuracy of the measurements of ZIMLAT 
and ESASDT of below one arcsecond. That the further 
observations in fact belong to the same object is validated 
via an orbit determination with both the observations used 
in the original sparse data orbit determination and the 
observations, which they were compared to. An orbit deter- 
mination with a root-mean-square of below two arcseconds 
is a reliable tool to associate observations of this accuracy 
of the same o bject to each other , as shown with cluster 
observations in iMusci et al] (l2005l ). 

2.2 Results 

In Fig. [1] the angular distance between predicted and 
observed position are displayed as a function of the time 
interval between the first and the last observation, which 
were used in orbit determination. Displayed are the mean 
values and the standard deviations of the angular distances 
of the single orbits. The mean value and standard devia- 
tions are determined with the single angular distances of 
predicted position to observed ones, all within 50 days since 
orbit determination. 

Figure [T] shows that the angular distances are in gen- 
eral very small. The vast majority of the determined orbits 
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to be interpreted as the effective area-to-mass ratio scaled 
with C/2 — 1; the AMR values of other sources may be 
scaled with a different factor. It is assumed that the AMR is 
constant over the orbital fit interval. A default value of 0.02 
m^kg~^ is selected, which corresponds to an AMR value 
of a standard GPS satellite, in case the AMR parameter 
is not estimated but kept fixed in the orbit determina- 
tion. For HAMR objects always an AMR value is estimated. 

The shadow paths of the orbit are modeled, under 
the assumption of a spherical earth on a mean circular 
orbit; the boundary between sunlit and eclipsed part is 
assumed to be cylindrical, no distinction between penumbra 
and umbra is made, earth atmosphere is neglected. 

For a long term investigation of the orbits and the 
AMR values, different comparable orbits have to be deter- 
mined. Only sparse observations are available, which are 
unequally spaced in time. A normalized setup is developed, 
tested with two low AMR objects and two of the HAMR 
objects of the AIUB catalog and applied for the creation 
of comparable orbits for the investigation of the HAMR 
objects. 

2 NORMALIZED SPARSE DATA SETUP 
2.1 The Method 

Four representative GEO objects from the internal cata- 
logue of the AIUB were chosen, they have been followed over 
longer time periods and are not listed in the USSTRAT- 
COM catalogue. Those objects are clearly space debris, 
since no maneuvers could be detected in the data. The 
AIUB did not have information what those objects actually 
were before becoming debris. From the apparent magnitude 
it can be concluded that those are all fragmentation pieces. 
They represent typical objects found in GEO surveys. Their 
properties are listed in Tab. [T] 

Two of the objects have low area to mass ratios, two objects 
qualify as HAMR objects with an AMR value larger than 
Im^ /kg. The optical angle-only observations are obtained 
with ZIMLAT (Zimmerwald, Switzerland), and ESASDT 
(Tenerife, Spain), supplemented by some observations of the 
ISON network provided by the Keldysh Institute of Applied 
Mathematics, Moscow, Russia. The latter observations were 
obtained from different sites of the ISON network, in these 
particular cases, all located in Eastern Europe. 

All orbits were determined from two observation sets 
only, using a priori orbital elements. A maximum of eight 
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Figure 1. Angular distances as a function of the time interval between the first and the last observation of the fit 
interval of orbit determination for object (a) E03174A, (b) E06321D, (c) E06327E and (d) E08241A. 
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Figure 2. Root mean square of orbit determination as a function of the arc length of observations for object (a) 
E03174A, (b) E06321D, (c) E06327E and (d) E08241A. 
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Figure 3. Angular distances as a function of the number of observations used for orbit determination for object (a) 
E03174A, (b) E06321D, (c) E06327E and (d) E08241A. 
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Figure 4. Angular distances as a function of anomaly distribution factor for object (a) E03174A, (b) E06321D, (c) 
E06327E and (d) E08241A. 
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Figure 5. Angular distance as a function of time interval covered by observations used for orbit determination for 
object (a) E03174A, (b) E06321D, (c) E06327E and (d) E08241A. 
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Figure 6. Time interval covered by the observations within the sets as a function of the number of observations used 
for orbit determination for object (a) E03174A, (b) E06321D, (c) E06327E and (d) E08241A. 



even produce distances smaller than 0.6 degrees. Except 
for the first object, each object also shows some outliers, 
with larger angular distances. These larger distances also 
tend to show larger standard deviations. The value of the 
angular distances seems to be, at least in this setup, quite 
independent of how large the difference between the first 
and the last observations of the fit interval is. Moreover, 
Figure [1] also shows that there is no significant difference in 
using observations only from one observation site for orbit 
determination or using observations from different sites. 
It could not be shown that the latter approach is more 
advantageous for orbit determination. Different observation 
sites still have advantages in terms of availability, weather 
conditions, which results in a larger amount of observations, 
which are available. In Fig. (2] the root mean square of the 
orbit determinations is shown, which were used for the 
prediction, as a function of the angular distance. No trend 
is visible, all orbits, which were determined had a small 
root mean square of below three arcseconds. 

In Fig. [3l the angular distances are displayed as a 
function of the actual number of single observations that 
entered orbit determination. It can be seen that no strong 
correlation is visible between the actual number of observa- 
tions used and the value for the distances. 

To find a measure for the true anomaly distribution, 
an anomaly distribution measure /ano was defined: It would 
be ideal to distribute all n observations equally spaced with 
an angle of 27r/n between each observation. The deviation 
from this ideal distribution is determined and normalized 
with the number of observations. The smaller /ano, the 



better distributed are the observations in anomaly. 



(— - (a^+i - ai)y + (— - (ai + 27r - an)) 

i—l 

(2) 

where as n is the number of observations and ai with 
i = l,..,n are the anomalies of the single observations, 
in ascending anomaly order. The angular distances as 
a function of /ano are displayed in Fig.|4l There is no 
clear correlation between the /ano and the distances, as 
it is expected for objects with small eccentricities. Object 
E06327E, with the highest eccentricity of e=0.06, has the 
strongest correlation with /ano- 

The crucial factor however, seems to be the time in- 
terval covered by observations itself within the sets. In 
Fig \5\ the angular distances are displayed as a function of 
the time interval covered within the two sets used in the 
beginning and the end of the fit interval, without the time 
gap in between the two sets. A strong correlation is visible. 
Fig. [6] shows that there is no strong correlation between 
the number of used observations and the time interval 
covered within the sets. For example for the ESASDT 
observation strategy, primarily densely spaced observations 
are available. 

The investigation of the data displayed in Fig. [6] showed 
that a coverage of at least 1.2 hours for both sets together 
seems to be necessary, in order to gain an orbit which allows 
to safely re-detect the investigated objects in more than 
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Figure 7. Inclination as a function of time for orbits of the object (a)E08241A, (b) E06321D, (c) E07194A, (d) E07308B, (e) E06293A. 
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Figure 8. Eccentricity as a function of time for orbits of the object (a) E08241A, (b) E06321D, (c) E07194A, (d) E07308B, (e) E06293A. 



90 percent of all cases with a field of view of one square 
degree, that is to have an accuracy of below 0.5 degrees. 



ranges between 10 and 120 days. As shown in the previous 
section the comparability of the orbits do not seem to be 
dependent on these ranges. 



3 INVESTIGATION OF HAMR OBJECTS IN 
SPARSE DATA SETUP 

The dynamical properties of HAMR objects were studied in 
the normalized sparse data setup established in the previous 
section. Orbits are determined with two observation sets 
only. The sets consist of four to eight observations each. The 
observations are required to span at least a time interval 
of 1.2 hours within the sets and need to be well spread 
over the anomaly for the objects in orbits with a high 
eccentricity. The total fit interval for orbit determination 



The orbits were first determined with observations from one 
observation site only, then with observations from different 
sites in the setup mentioned above. The observations used 
in this investigation stem from the ESASDT, ZIMLAT, and 
from several telescopes of the ISON network. 
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Figure 9. AMR as a function of time for orbits of the object (a) E08241A, (b) E06321D, (c) E07194A, (d) E07308B, (e) E06293A. 
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Figure 10. Relative variation of AMR value as a function of the absolute AMR value of (a) 47 HAMR and (b) LAMR objects. 



3.1 Selected Objects 

Five objects were selected for a detailed investigation. All 
objects were discovered and first detected by the AIUB and 
are not hsted in the USSTRATCOM catalogue. All objects 
are faint debris objects. They were tracked successfully over 
several years, and no maneuvers were detected. A set of 
osculating orbital elements and an average value for the 
apparent magnitudes are listed in Tab.O The two objects 
with the lowest AMR values, E08241A and E06321D, which 
were used in the investigation of the sparse data orbit 
determination, are used here again. 



3.2 Evolution of Orbital Elements 

The evolution of the orbital elements over time is in- 
spected in a first step. Figure[7l shows the development 
of the inclination and its errors in inclination, of the 
five objects. The error bars are too small, to be visible 
in the plot in most cases. The inclination values of the 



Table 2. Investigated HAMR objects: Internal name, 
epoch (MJD), eccentricity, inclination (deg), semi- major 
axis (km), area to mass ratio {m? /kg) and apparent mag- 
nitude (mag). 



NAME 


Epoch 


a 


e 


i 


AMR 


Mag 


E08241A 


55213.0 


41600 


0.041 


13.26 


1.24 


16.1 


E06321D 


55275.9 


41400 


0.035 


7.00 


2.29 


15.3 


E07194A 


54877.0 


40900 


0.005 


7.31 


3.37 


16.8 


E07308B 


54416.0 


35600 


0.264 


7.63 


8.83 


15.8 


E06293A 


54951.0 


40200 


0.245 


11.06 


15.41 


16.8 



different orbits are closely aligned to each other and mark 
a consistent evolution, only in the case of object E08241A 
in Fig. [71 a wider spread in the inclination values can be 
observed. The orbits determined with observations from 
the different observation sites produce almost identical 
results. For object E07308B and E06293A, which have the 
highest AMR values, the inclination seems not to follow 
a steady increase over time, but some smaller periodic 
substructure seems to be superimposed. These may very 
well be the perturbations with a period of one nodal 
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Figure 11. Angular distances of predicted orbits on the celestial sphere as a function of AMR for orbits of the object (a) E08241A, (b) 
E06321D, (c) E07194A, (d) E07308B, (e) E06293A. 
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Figure 12. Error of the AMR value as a function of AMR as estimated in orbits of the object (a) E08241A, (b) E06321D, (c) E07194A, 
(d) E07308B, (e) E06293A. 



year, w hich are well known f or objects with hig h AMR , 
seee.g.. [Liou Weaverl (|2005l ). ISchildknecht et a/. 1 (|2005bh . 

Figure [8] shows the evolution of the eccentricity values 
and its errors estimated in orbit determination for the 
different objects. Periodic variations can be observed for all 
objects. The different orbits with observations from one site 
only or from different sites result in the same eccentricities. 



3.3 Evolution of Area to Mass Ratio Value 

Figure [9] shows the AMR values as a function of time for 
the objects listed in Tab. [2] In all cases, the values for the 
AMR do not show clear and obvious common trends, see 
Fig.[Zland[8l 

For object E08241A, the AMR values vary around a 
mean value of 1.4m^kg~^ with no obvious trend or periodic 
signal, see Fig. [9^. 



For object E06321D (see Fig.Eb), the AMR value seems to 
vary periodically with a period of about one year around 
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Figure 13. Absolute values and standard deviations of the angular distances as a function of the error of the AMR value as found in 
orbit determination of the object (a) E08241A, (b) E06321D, (c) E07194A, (d) E07308B, (e) E06293A. 



a value of 2.5m^kg ^, but also values of 2.35 m^kg ^ 
and 2.65m^kg;~^ o ccur. Similar results were obtained by 
iMusci et a/~l (|2010l ). for the same object, in different orbit 
determination setups. The AMR value of object E07194A 
(see Fig.[9t) varies around 3.5m^kg~"^, but in the orbits 
determined with combined observations from all the sites, 
so-called outliers of 4.5m^kg~"^ and 2.3m^kg~"^ occur as 
well. These have, however large error values. 

Object E07308B (see Fig.EJi) seems to generally in- 
crease its AMR value over time from a value of 8.5m^kg~^ 
up to 9.0m^kg~^. But single orbits also show AMR values 
of i.e. lOm^kg"^ 

Figure [9^ shows that object E06293A, which is the 
object with the largest AMR value regarded here, has 
significant data gaps. A general trend of the AMR value in 
time, increasing from 15.5 m^kg""*^ to 16.5 m^kg""*^ cannot 
be excluded. But one orbit determined with ESASDT data 
also shows a value of 18.2 m^kg~^, with a small formal error. 

No general correlation between the AMR value itself 
and the variations of the AMR value could be determined, 
no general trend is visible. A study on the variation 
of AMR values wa s conducted by ISchildkecht et al~\ 
(|to be published 2Qllh . The variations of the AMR values 
of 47 HAMR objects were investigated and compared to 
the AMR variations of orbits of 40 low AMR (LAMR) 
value objects. No normalized or sparse data setup orbit 
determination setup was chosen. The AMR values in that 
analysis were determined in the standard orbit determina- 
tion procedure at the AIUB, with fit arcs as long as possible 
for a successful, that is defined as leading to a small rms 
error, orbit determination. The results are illustrated in 
Fig. 1101 No general trend in the AMR variations could 
be determined for either HAMR or LAMR objects. The 
relative variations of the AMR values of the LAMR objects 



were larger, than the AMR variations of the HAMR values. 
The AMR variations of the LAMR objects were of the 
order of several 100 percent. 

All orbits were predicted and compared to additional 
observations of the same object, which were not used for 
orbit determination. The additional observations were all 
checked via orbit determination, to ensure that they belong 
to the same object. FigurefTTl shows the angular distances 
between the predicted ephemeris and observations. The 
values are averaged over all distances 50 days after orbit 
determination and their standard deviations serve as error 
bars. 

FigurefTTk shows that for object E08241A, one orbit 
produces the largest distances of one degree. This orbit 
does not show up prominently in the orbital parameter 
plots (see Fig.fSK and[7K) or AMR value plots (see Fig. [9^). 
The orbit with ZIMLAT data, which produced the outlier 
AMR value of 0.82m^kg~^, does not show up prominently 
in the distance plot fFig.lllb). 

The mean value of all angular distances of object E06321D 
are well below 0.2 degrees, but four orbits show large 
standard deviations in the angular distance, as Fig.fTTb 
shows. All of them have been determined with combined 
observations from ZIMLAT, ESASDT, and ISON obser- 
vations. Their AMR values are 2.36 m'^kg~^, 2.50m^kg"\ 
2.57m2kg"\ and 2.66 m^kg"^ The orbits with the AMR 
value of 2.36m^kg~^ does show up also in a group of 
outlier AMR values, which do not seem to follow the 
periodic variation in the evolution of the AMR values. The 
other orbits, with large standard variations in the angular 
distance do not show up prominently (Fig.[9l3). Those orbits 
with the largest standard variation in angular distance do 
not show the largest error in the AMR values either, as 
Fig. Us] shows. 
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FigurefTTt shows for object E07194A three angular 
distances with large standard deviations. The orbits were 
determined with observations from all sites. They have 
AMR values of 2.12m2kg"\ 2.21m2kg"\ and 4.46 m^kg"^ 
Those are the smallest and largest AMR values in the deter- 
mined orbits for E07194A. These three values do also show 
up as outliers in Fig.[9]:. For objects E07308B and E06293A, 
the angular distances with a large standard variation (see 
Fig. Hid and e), do not show significant outlier AMR values 
in Fig.EJi and e. For object E07308B, the orbit with an 
AMR value of 10.15 m^kg~^ shows the largest mean value 
in the angular distance of almost 0.7 degrees but has a small 
standard deviation in this distance (Fig.[TTJl). This value 
is significantly different compared to the other determined 
AMR values, see Fig.[9li. 

The dependency of the AMR value on the error of 
the AMR, as it was found in orbit determination, is 
investigated in the final step. No clear correlation could 
be determined between the AMR value and its rms value 
(Fig.[l2)). 

Figure[T3] shows the angular distance distances on the 
celestial sphere as a function of the error of the AMR value. 
As expected, for none of the objects a clear correlation 
between the error of the AMR value and the absolute value 
of the distances or the standard deviation of the distances 
could be determined. 

All investigated objects show variations in the AMR 
value, but not a common characteristic in these variations. 
It has to be noted that the result may be affected by the 
relatively simple shadowing m odel that was used i n orb it 
determination; as it wa s shown [Pardini Anselrnol (| 20081 ). 
IValk Lemaitrd (|2008l ). shadowing effects have a signifi- 
cant influence on the long term evolution of orbits of HAMR 
objects. Investigation of simulated orbit s with numerical 
and se mi- analytica l methods, e.g. by IValk Lemaitrel 
(I2OO8I ). IValk et al I (I2OO8I ) also showed the existence of 
irregular chaotic orbits and the significant influence of 
secondary resonances on the orbits of HAMR objects. 
However, these simulations did assume a constant AMR 
value. Complex attitude motion, irregular shapes, and/or 
deformation of the actual objects, could lead to an actual 
change in the AMR value itself over time, which may not 
be averaged out over the fit interval of orbit determination. 



4 CONCLUSIONS 

A sparse data setup was established to create compa- 
rable orbits over longer time intervals. Orbits with two 
data sets only do produce small differences between the 
propagated ephemerides and further observations, as long 
as 1.2 hours are covered within the sets. Other factors, 
such as that the observations stem from different sites or 
the time interval between the sets, are found to be negligible. 

The orbits of high area-to-mass ratio (HAMR) objects 
were analyzed in this setup. The AMR value, that is 
the scaling factor of the direct radiation pressure (DRP) 



parameter, varies over time. The order of magnitude of the 
variation of the area-to-mass ratio (AMR) value was not 
correlated with the order of magnitude of its error. 

The variation of the AMR is not averaged out in the 
fit interval of orbit determination. In the evolution of the 
AMR value over time, no common characteristic could 
be determined for different HAMR objects. Further work 
on the orbits of HAMR objects is needed, to improve the 
radiation pressure model, to determine possible attitude 
motion or deformations and to understand also resonance 
effects and the existence of chaotic regions. 
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